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1.  STATEMENT  OF  PROBLEM 

Mathematical  models  for  describing  the  failure  of  materials  at  high  rates  of 
loading  arc  inadcouate  for  computing  the  response  of  structures  subjected  to  such 
loading.  For  example,  improved  models  arc  required  for  describing  the  development 
of  bands  of  localized  shearing  deformation  ("shear  bands")  which  form  in  thick  steel 
plates  when  these  plates  arc  impacted  by  projectiles.  Improved  models  are  required 
also  for  describing  the  fracture  of  metals  under  impact  loading  conditions.  Such 
models  should  be  based  on  a  firm  understanding  of  the  failure  mechanisms  in  order 
to  determine  which  material  properties  must  be  measured  to  predict  the  response  of 
the  structure  over  its  range  of  loading  conditions.  Furthermore,  a  mechanistic 
understanding  is  required  in  order  to  guide  the  development  of  new  materials  and 
microstructures  with  improved  properties  at  high  loading  rates. 

Under  this  research  agreement,  this  overall  problem  has  been  addressed  by  five 
faculty  investigators  who  have  investigated  critical  conditions  for  dynamic  fracture 
initiation  and  for  the  formation  of  shear  bands.  These  investigations  have  been 
supplemented  by  investigations  of  the  plastic  response  of  metals  at  the  high  strain 
rates  that  occur  in  shear  bands  and  in  crack-tip  zones.  Analytical  and  computational 
approaches  have  been  developed  to  interpret  the  various  experiments  and  to  provide 
understanding  for  the  development  of  mechanistic  models  for  computing  the  failure 
of  structures  subjected  to  high-rate  loading. 

Research  results  obtained  under  the  agreement  arc  described  in  the  following 
section.  References  without  superscripts  refer  to  publications,  reports  and  theses  that 
were  supported  by  the  current  grant;  these  publications  are  listed  in  Section  2.4. 
References  not  supported  by  the  grant  arc  denoted  by  superscripts  R,  c.g.  ( l ‘>87 
and  listed  in  Section  2.5. 
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2.  SUMMARY  OF  IMPORTANT  RESULTS 


2.1  Shear  Rands 


2.1.1  Experiments 

A  torsional  Kolskv  bar  (split-Hopkinson  bar)  has  been  used  to  deform 

tubular  specimens  of  four  steels  (cold  rolled  A1SI  1018,  hot  rolled  AISI  1010,  4340 

VAR,  and  HY-100)  at  strain  rates  of  103  s'1  [HARTLEY,  DUFFY  and  HAWLEY, 
1987],  Shear  bands  have  been  observed  to  form  in  all  of  these  steels.  The 

temperature  of  the  material  in  the  bands  has  been  measured  by  monitoring  the 
infrared  radiation  emitted  at  the  metal  surface.  A  linear  array  of  ten 

indium-antimonide  detectors  has  been  used  to  monitor  the  temperature  history  at  ten 
neighboring  points  across  the  gage  section  of  the  specimen.  For  the  relatively  wide 
(150/im)  shear  bands  obtained  for  the  AISI  1018  and  AISI  1020  steels  the  maximum 
temperature  rise  in  the  band  is  approximately  450°C.  This  value  agrees  well  with  the 
temperature  rise  estimated  for  adiabatic  deformation  of  the  material  to  the 

shear-band  strains  observed  through  the  distortion  of  a  grid  copied  onto  the  specimen. 
The  shear  strains  at  localization  arc  approximately  15%  for  the  cold  rolled  steel  and 
100%  for  the  hot  rolled  steel.  This  large  difference  in  the  critical  nominal  strains  at 
localization  is  consistent  with  predictions  that  account  for  the  much  greater  strain 
hardening  observed  for  the  hot  rolled  steel. 


In  order  to  measure  the  temperature  increases  in  the  relatively  narrow  bands 
(-  10(zm)  obtained  for  the  4340  VAR  and  HY-100  steels,  the  infrared  detection 
system  was  improved  to  minimize  aberrations.  With  this  improvement,  and  with  spot 
sizes  of  35^01,  temperature  increases  up  to  600°C  have  been  recorded.  From  the 
strains  estimated  by  observing  the  distorted  grids,  the  estimated  temperatures  in  the 
shear  bands  arc  expected  to  be  larger  than  the  recorded  values.  This  discrepancy  is 
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attributed  to  the  lack  of  sufficiently  fine  resolution  in  the  infrared  detection  system 
to  allow  the  accurate  measurement  of  temperature  over  such  narrow  bands. 
Availability  of  these  temperature  measurements,  combined  with  new  high-speed 
photographs  of  the  developing  shear  band,  has  enabled  us  to  establish  that  the 
formation  of  the  intensely  localized  shear  band  coincides  with  the  sharp  decrease  in 
the  torque  transmitted  through  the  specimen.  The  critical  strain  at  which  the  drop 
in  load  carrying  capacity  occurs  varies  from  specimen  to  specimen  in  a  manner  that 
appears  to  be  consistent  with  the  predicted  effects  of  geometric  imperfections  on 
shear  strain  localization  in  simple  shear  [MOL1NARI  and  CLIFTON,  1987], 

2.1.2  Theory 

Shear  band  formation  during  simple  shearing  of  thermoviscoplastic  materials 
has  been  investigated  by  means  of  either  closed  form  or  essentially  closed  form 
solutions  of  a  simplified  model  problem.  In  the  model  problem,  inertia,  heat 
conduction,  and  clastic  deformations  aTe  neglected.  Full  numerical  solutions  have 
shown  that  these  effects  arc  essentially  negligible  for  the  strain  rates,  specimen  sizes, 
and  materials  used  in  the  torsional  Kolsky  bar  experiments.  Previously,  MOLINARI 
and  CLIFTON  ( 1 983 )R  had  shown  that,  for  the  model  problem,  and  for  several 
elementary  constitutive  models,  simple  analytical  expressions  could  be  obtained  for  the 
critical  conditions  of  strain  rate  sensitivity,  thermal  softening  and  strain  hardening 
for  which  shear  localization  would  occur  for  sufficiently  large  strains.  These  results 
have  been  extended  to  account  for  geometric  and/or  material  imperfections.  For 
given  imperfections  the  critical  nominal  shear  strain  at  localization  has  been  obtained 
either  in  closed  form  or  in  the  form  of  integrals  that  are  relatively  easy  to  evaluate 
[MOLINARI  and  CLIFTON.  1987],  Wall  thickness  variations  of  specimens  used  in  the 
torsional  Kolsky  bar  experiments  have  been  measured  and  used  in  the  prediction  of 
the  relationship  between  the  nominal  shear  stress  and  the  nominal  shear  strain  during 


shear  band  lormation.  These  predicted  stress-strain  curves"  show  good  qualitative 
agreement  with  those  measured  in  experiments;  quantitative  agreement  is  sufficiently 
close  to  indicate  that  accurate  measurement  of  geometrical  imperfections  of  specimens 
is  important  in  making  detailed  comparisons  between  theory  and  experiment  in  shear 
band  experiments. 


Investigation  of  the  effects  of  strain  hardening,  strain-rate  sensitivity,  thermal 
softening,  heat  conduction  and  the  imposed  strain  rate  on  the  shear  localization  process 
has  been  extended  to  the  ease  of  plane  strain  compression  [LEMONDS  and 
NEEDLEMAN,  1986a],  The  deformation,  stress  and  temperature  fields  arc  computed  in 
an  infinite  solid  which  contains  a  periodic  rectangular  array  of  inhomogcneitics.  The 
inhomogcncitics  give  rise  to  non-uniform  deformation  fields  which  may  localize  in  the 
form  of  a  shear  band.  Because  these  fields  arc  periodic,  only  a  rectangular  region 
surrounding  a  single  inhomogcncity  needs  to  be  considered.  Full  two-dimensional,  finite 
element  computations  arc  performed  within  the  context  of  a  viscoplasticity  theory  which, 
in  the  rate  independent  limit,  corresponds  to  flow  theory  with  combined  isotropic  and 
kinematic  hardening.  Full  account  is  taken  of  finite  strain  and  rotation  effects,  but 
attention  is  confined  to  quasi-static  loading.  The  predicted  response  depends  significantly 
on  the  multi-axial  hardening  characterization  of  the  solid.  The  initiation  and 
propagation  of  shear  bands  is  much  more  pronounced  for  the  limiting  case  of  kinematic 
hardening  than  for  the  other  limiting  ease  of  isotropic  hardening.  This  difference  in 
the  susceptibility  of  various  model  materials  to  shear  band  formation  is  viewed  as  an 
indication  that  constitutive  models  involving  flow  potentials  with  large  curvature  near 


the  current  stress  state  (as  for  corner  theories  in  rate  independent  plasticity)  are  models 
for  which  predictions  of  shear  strain  localization  arc  more  likely. 


In  order  to  develop  a  better  understanding  of  the  effects  of  thermal 
softening  and  heat  conduction  in  shear  localization,  these  effects  arc  included  in  the 
analysis  of  localization  in  an  infinite  band  [LEMONDS  and  NEEDLEMAN,  1986b], 
Full  account  is  taken  of  finite  geometry  changes,  but  inertial  effects  arc  neglected. 
An  energy  balance  is  written  between  homogeneously  deforming  bands  in  a  manner 
that  models  conditions  in  the  above  finite  clement  study  of  shear  band  development. 
The  infinite  band  analysis  requires  specification  of  an  imperfection  amplitude  and  a 
length  scale  over  which  heat  conduction  effects  are  significant.  These  are  chosen  to 
match  results  of  the  finite  element  analysis.  The  predictions  of  the  band  analysis 
and  the  full  finite  element  computations  arc  then  compared  for  a  wide  range  of 
material  properties  for  both  isotropic  and  kinematic  hardening  characterizations  of 
the  flow  potential  surface.  The  predicted  dependence  of  the  onset  of  localization  cn 
material  properties,  such  as  strain  hardening  and  strain-rate  sensitivity,  is  the  same 
for  both  types  of  analyses. 

The  effects  of  inertia  on  shear  band  development  in  plain  strain  compression 
ha\e  been  investigated  for  an  elastic-!  iscoplastic  material  with  softening  introduced 
through  a  hardness  function  that  exhibits  a  local  maximum  [NEEDLEMAN,  1987], 
An  initial  inhomogcncity  is  introduced  to  cause  the  deformation  to  be  non-uniform. 

Regardless  of  whether  the  material  is  hardening  or  sol'tening,  plastic  strain 

development  involves  the  evolution  of  finger-like  contours  emanating  from  the 
inhomogcncity  at  45°  to  the  compression  axis.  Once  a  given  strain  contour  crosses 
the  specimen,  it  fans  out  about  its  initial  direction  of  propagation.  For  a  softening 
solid,  this  fanning  out  ceases  for  some  .train  level  greater  than  the  strain  at  the 
hardness  maximum  and  further  aiuimr.g  takes  place  in  an  ever  narrowing  band. 

Many  of  the  qualitative  feature  of  .car  land  development  under  dvn.miic  loading 
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Dynamic  Fracture 


2.2.1  Notched  Bar  Experiments  and  Interpretation 

An  investigation  has  been  conducted  into  the  effects  of  microstructurc, 

loading  rate,  and  temperature  on  the  initiation  of  plane  strain  fracture  of  a  plain 

rate,  carbon  A1SI  1020  steel  [COUQUE,  ASARO,  DUFFY,  and  LEE,  1987],  Ferrite 
grain  sizes  and  prior  austenite  grain  sizes  were  varied  by  changing  the  austenitizing 
temperature  and  the  cooling  rates  during  normalization.  Three  different 
microstructures  were  investigated.  Grain  size  differences  in  these  three  microstructures 
were  designed  to  provide  independent  variation  of  each  of  the  grain  sizes  while 
varying  these  sizes  by  a  factor  of  2  or  more.  Fracture  toughness  tests  were 

conducted  using  pre-crackcd  notched  round  bars  loaded  in  tension.  The  tests  were 

u 

performed  at  quasi-static  loading  rates  Kj  =  1  MPa  m^s  ,  and  at  dynamic  loading 

rates  Kj  =  2  x  106  MPa  m^s'1.  Testing  temperatures  covered  the  range  from  -150°C  to 
150°C  corresponding  to  fracture  initiation  modes  varying  from  transgranular  cleavage 
to  fully  ductile  fracture.  In  addition,  Charpy  impact  tests  were  conducted  over  the 
same  range  of  temperatures.  The  plastic  shearing  response  of  the  steels,  over  this 
range  of  temperatures,  was  obtained  through  torsion  tests  at  quasi-static  shear  strain 
rates  of  7  =  5.0  x  10'4s'1  and  at  dynamic  shear  strain  rates  of  j  =  1.5  x  lOV1. 

The  measured  dependence  of  fracture  toughness  on  microstructurc,  loading 
rate  and  temperature  is  shown  in  Fig.  1.  The  microstructurc  T!  has  the  smallest 
grain  sizes  and  the  microstructurc  T7  has  the  largest  grain  sizes;  T9  has  the  prior 
austenite  grain  size  of  T1  and  a  ferrite  grain  size  comparable  to  that  of  T7.  For  all 
three  microstructuies  the  temperature  dependence  of  the  dynamic  I  racturc  toughness 

Kjr)  consists  of  a  "lower  shelf"  at  temperatures  below  room  temperature  iRD,  an 
"upper  shelf"  at  temperatures  above  100"C,  and  a  transition  over  a  relatively  narrow 
range  of  temperatures  between  RT  anil  IOOr'C.  The  temperature  dependence  ol  the 
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quasi-static  fracture  toughness  K.[c  is  more  complicated.  Although  for  each  of  the 
three  microstructurcs  Kk  increases  from  a  low  value  at  150'’C  to  an  apparent  "upper 
shelf”  at  temperatures  above  50°C.  the  increase  is  not  monotonic  and  the  "two  shelf" 
structure  is  not  clearly  delineated.  For  all  microstructurcs  and  loading  rates,  fracture 
initiation  at  -150°C  occurs  by  cleavage  and  fracture  initiation  at  +150°C  occurs  by 
ductile  rupture  processes.  The  fracture  toughness  at  -150°C  is  essentially  the  same 
for  all  microstructurcs  and  loading  rates.  At  +150°C  the  fracture  toughness  increases 
with  increasing  loading  rate  and  varies  significantly  with  variations  in  the 
microstructurcs.  These  variations  in  fracture  toughness  for  fully  ductile  fracture 

processes  follow  the  variations  in  flow  stress  with  loading  rate  and  microstructurc. 

In  order  to  get  a  better  understanding  of  the  relationship  between  the 
fracture  toughness  and  the  micromechanical  processes  involved  in  fracture  through 
cleavage  and  through  void  growth,  the  fractured  specimens  were  examined  using 

transmission  and  scanning  electron  microscopy.  For  cleaved  surfaces  the  size 
distributions  of  cleavage  facets  were  measured.  Dimple  size  distributions  were 
obtained  for  surfaces  that  had  separated  by  void  growth  and  coalescence.  These 

measurements  were  used  with  micromcchanical  fracture  models  in  an  attempt  to 
correlate  the  measured  fracture  toughnesses  with  parameters  that  characterize  the 
microscopic  failure  processes.  For  pure  cleavage,  fracture  was  assumed  to  occur  when 
the  maximum  principal  stress  exceeds  a  critical  value  over  a  critical  distance  from 
the  crack  tip.  The  critical  distance  was  taken  to  be  two  ferrite  grain  diameters  --  a 
distance  that  correlates  well  with  cleavage  facet  sizes.  For  fully  ductile  rupture, 
fracture  was  assumed  to  occur  when  an  equivalent  plastic  strain  exceeds  a  critical 

value  over  a  critical  distance  from  the  crack  tip.  The  critical  strain,  obtained  from 
independent  tensile  tests  on  notched  bars,  was  taken  to  be  the  strain  at  which  void 
nuclcation  occurs.  The  critical  distance  was  taken  to  be  twice  the  spacing  id 


fracture  initiation  sites,  e.g.  MnS  inclusions,  pearlite  grains,  and  grain  boundary 
carbide  plates.  For  both  crack  propagation  models  the  dependence  of  stress  and 
strain  on  radial  distance  from  the  crack  tip  was  taken  from  the  finite  element 
solutions  of  McMEEKING,  (I977)R.  Comparison  of  the  calculated  and  measured 
fracture  toughness  values  arc  shown  in  Table  1.  The  agreement  is  quite  good,  even 
on  the  non-monotonic  increase  in  K.Ic  with  test  temperature  found  for  the 
microstructurc  Tl.  This  agreement  suggests  that  quantitative  micromodeling  using  a 
simple  fracture  model  is  a  promising  means  of  correlating  both  quasi-static  and 
dynamic  fracture  toughness  with  microstructures. 
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Tabic  I  Calculated  and  measured  fracture  toughness  values. 


Tl 

T7 

Loading 

Condition 

Testing 

Temperature 

Calculated 

KIc’KId 

Measured 

KIc,KId 

Calculated 

KIc*KId 

Measured 

KIc’KId 

[°C] 

[Mpa  /  m] 

[Mpa  / mj 

[Mpa  / mj 

[Mpa  v  'mj 

Static 

-  150 

30  (C) 

20  (C) 

39  (C) 

27  (C) 

-  60 

89  (F) 

91  (F) 

57  (Cj 

57  1C) 

23 

73  (F) 

70  (F) 

7S  iF) 

94  (FI 

100 

54  (F) 

101  (F) 

9S  (F) 

i 

111 iF) 

Dynamic 

-  150 

25  (C) 

21  (C) 

3S  «C) 

22  (O 

-  60 

40  (C) 

31  (C) 

44  (C) 

3S  iC) 

23 

5S  (C) 

54  (C) 

57  (C) 

51  O 

100 

120  IF) 

14S  (F) 

1 60  id 

135  d 

2.1.2  Dynamic  Fracture  in  Plate  Impact  Experiments 

A  new  plate  impact  experiment,  originated  under  the  previous  grant,  has  been 
developed  for  the  study  of  dynamic  fracture  processes  that  occur  on  the  sub-microsecond 
time  scale  [R AVICHANDR AN  and  CLIFTON,  1987],  The  experiment,  designed  to  make 
the  measurements  amenable  to  interpretation  within  the  framework  of  fracture 
mechanics,  corresponds  to  the  plane-strain  loading  of  a  semi-infinite  crack  by  a  step 
tensile  pulse  at  normal  incidence.  To  obtain  these  conditions  a  disc  containing  a 
mid-plane,  pre-fatigued,  edge  crack  that  has  been  propagated  halfway  across  the  diameter 
is  impacted  by  a  thin  flyer  plate  of  the  same  material.  A  compressive  pulse  propagates 
through  the  specimen  and  reflects  from  the  rear  surface  as  a  square  tensile  pulse  with  a 
duration  of  l^s.  This  plane,  tensile  pulse  causes  initiation  and  propagation  of  the 
fatigue  crack.  The  crack  advance  under  the  known  pulse  is  measured  by  using  a 
focussed  ultrasonic  transducer  to  locate  the  crack  front  before  and  after  the  plate  impact 
experiment.  Various  models  for  the  advance  of  the  crack  arc  examined  by  prescribing 
the  crack  tip  motion  corresponding  to  these  models  and  comparing  the  predicted  motion 
at  a  point  on  the  rear  surface  of  the  specimen  with  the  motion  monitored  during  the 
experiment  using  a  laser  velocity  interferometer.  All  comparisons  arc  made  before  the 
experimental  results  can  be  affected  by  the  arrival  of  unloading  waves  from  the 
periphery  of  the  specimen. 

Experiments  have  been  conducted  on  a  hardened  4340  VAR  steel  (Rc  =  52) 
at  room  temperature  and  at  temperatures  of  approximately  -100°C.  Crack  advance 
increases  with  increasing  impact  velocity  and  decreasing  temperature.  The  crack 
paths  arc  quite  straight,  with  essentially  no  branching;  significant  turning  of  the 
crack  towards  the  rear  surface  of  the  specimen  occurs  at  late  times.  Examination  of 
the  fracture  surfaces  by  means  of  scanning  electron  microscopy  indicates  that  at  room 
temperature  fracture  occurs  primarily  by  the  ductile  process  of  void  initiation,  growth 


nnd  coalescence  whereas  at  low  temperatures  fracture  occurs  by  a  mixed  process 
in\ Giving  cleavage  across  irregularly-shaped  regions  with  a  diameter  of  approximately 
lOjcm  combined  with  ductile  rupture  of  smaller  regions  between  the  cleavage  facets. 
The  dynamic  fracture  toughness,  represented  by  the  critical  value  of  the  stress 
intensity  factor  K,c  required  for  fracture  initiation,  is  inferred  from  the 
corresponding  elastodv namic  solution  and  the  measured  advance  of  the  crack.  Values 
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Figure  2  Correlation  between  measured  and  predicted  crack  advances  for 
plate  impact  experiments  on  4340  VAR  steel  (R.  =  52). 

Predictions  arc  shown  for  three  crack  propagation  models:  (i) 
constant  energy  release  rate.  G  =  Gc;  (ii)  constant  stress  intensity 
factor  K  =  KIc;  (iii)  constant  crack  velocity,  a/cR  =  0.4,  0.5,  0.6. 
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ot  K,c  =  62  MPam  and  K.,  =  32  MPam  arc  obtained,  respectively,  lor  tests  at  room 
temperature,  and  tests  at  -I00°C.  Essentially  the  same  values  arc  obtained  for  three 
different  crack  propagation  models:  constant  energy  release  rate,  constant  stress 
intensity  factor,  and  constant  crack  velocity  at  approximately  one-half  the  Rayleigh 
wave  speed  (see  Fig.  2);  for  all  models  the  critical  stress  intensity  factor  at  crack 
arrest  is  assumed  to  be  the  same  as  at  crack  initiation.  The  values  for  K,c  arc 
higher  than  interpolated  values  (K.Ic  =  46  and  KIc  =  21,  respectively)  obtained  from 
quasi-static  notched  bar  experiments  [CHI,  LEE  and  DUFFY,  1987]  on  a  steel  with 
the  same  chemical  composition.  However,  the  values  for  K.Jc  arc  lower  than  the 
values  (Klc  =  69  and  K,c  =  54,  respectively)  obtained  in  dynamic  notched  bar 
experiments  [CHI,  LEE»  and  DUFFY,  1987]  at  loading  rates  of  Kj  =  106  MPam^s'1. 
The  lower  values  of  K.Ic  obtained  at  the  extremely  high  loading  rates  (i.c.  K,  =  2  x 
I08  MPam'V1)  of  the  plate  impact  experiments  have  recently  been  attributed  [GODSE, 
RAVTCHANDRAN  and  CLIFTON,  1988]  to  the  presence  of  packets  of  upper  bainite 
which  fail  by  cleavage  at  load  levels  below  those  which  cause  failure  of  the  fully 
martensitic  microstructurc.  The  presence  of  upper  bainite  (up  to  15%)  results 
because  the  plate  impact  specimens  arc  cut  from  2.5-inch  diameter  bars  which  are  too 
large  to  allow  the  fast  quenching  at  the  center  that  is  required  to  obtain  a  fully 
martensitic  microstructurc. 

A  second  order  accurate,  finite  difference  method  has  been  developed  for  the 
numerical  simulation  of  these  experiments.  flic  loading  is  modeled  as  that  of  a  plane, 
square,  tensile  pulse  impinging  at  normal  incidence  on  a  semi-infinite  crack  in  a  plate 
bounded  by  two  free  surfaces.  Crack  advance  is  assumed  to  initiate  when  the  stress 
intensity  factor  reaches  the  critical  value  K.,.  obtained  from  correlations  with  the 
measured  total  crack  advance.  For  the  low  temperature  experiments  in  which  the 
crack-tip  plasticity  appears  to  be  sufficiently  limited  to  justify  the  use  ol  linear 
clastic  fracture  mechanics,  tire  crack  velocities  have  been  prescribed  according  to  the 
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three  different  crack  propagation  models  mentioned  in  the  previous  paragraph. 
Although,  as  noted  previously,  the  value  of  obtained  from  the  measured  total 

crack  advance  appears  to  be  relatively  insensitive  to  the  crack  propagation  model,  the 
predicted  motion  at  the  rear  surface  is  quite  sensitive  to  the  crack  propagation  model. 
Best  agreement  with  the  measured  motion  at  the  rear  surface,  during  the  motion  due  to 
the  initial  step  wave,  appears  to  be  obtained  for  the  constant  crack  velocity  model  (sec 
Fig.  3).  For  the  room  temperature  experiments,  viscoplastic  effects  arc  included  in 
the  computations.  Reasonable  agreement  in  predicted  and  measured  free  surface 
velocity-time  profiles  is  again  obtained  when  the  crack  velocity  is  taken  to  be  constant. 


-  Constant  Crack  Velocityi(a/  CR;  =  0.6) 

. .  Constant  Energy  Release  Rate 

Constant  Stress  Intensity 
- -  Stationary  Crack 

O  Experiment  i^Shot:  35-20) 
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Figure  3  Comparison  si'  measured.  f roe-surface,  \  eloeity -time  profile 
for  a  low  temperature  ( - 1 00"C )  plate  impact  experiment  with 
profiles  predicted  for  three  crack  propagation  models:  (i) 
constant  energy  release  rate:  (ii)  constant  stress  intensity 
factor:  (iii)  constant  crack  \clocttv.  Fhe  predicted 

velocitv-timc  profile  for  a  national  \  track  is  included  to 


(actor:  (iii)  constant 

velocity-time  profile  for 
show  the  .en.iti\it\  of 
motion  of  the  crack. 
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2.2.3  Dynamic  Crack  Propagation  Theory 

A  model  ol  the  high  strain  rate  crack  growth  process  'a  as  dev  eloped  in 
j  collaboration  with  Professor  J.W.  Hutchinson  of  liars  ard  I’nis  ersity.  based  ■  n  •  es  oral 

observations.  First,  estimates  of  material  strain  rates  near  the  tip  of  a  crack 
growing  in  an  clastic-plastic  material  under  common  conditions  suggest  that  values  in 
I  the  range  I04  -  10°  per  second  may  be  anticipated.  Fxperiments  on  iron  conducted 

at  such  rates,  under  more  or  less  homogeneous  states  of  deformation,  indicate  that  the 
high  strain  rate  response  may  be  divided  into  two  major  regimes.  For  strain  rates 

|  below  a  certain  level,  the  flow  stress  is  relatively  insensitive  to  changes  in  strain 

rate,  whereas  above  this  level  the  dependence  is  quite  strong.  1  he  problem  is  then 
studied  from  the  small  scale  yielding  point  of  view  ol'  nonlinear  fracture  mechanics. 

1  Crack  growth  is  assumed  to  occur  in  a  bod>  that  remains  elastic  at  points  remote 

from  the  crack  tip.  The  applied  loading  is  icflcctcd  in  a  remote  stress  intensity 
factor  or,  equivalently,  a  rate  of  energy  flux  into  the  crack  tip  region  per  unit  crack 
advance.  In  the  crack  tip  region,  the  potentially  large  stresses  are  relieved  to  some 
extent  through  viscoplastic  flow.  Thu  .  crack  growth  is  accompanied  by  a  /one  of 
active  plastic  deformation  and  a  wake  region  >,f  permanently  deformed  but  unloaded 
material  along  the  crack  flanks.  Within  the  inner  portion  of  the  active  plastic  zone, 
where  the  strain  rates  arc  such  that  the  response  is  in  the  high  rate  regime,  the 

stress  field  is  also  singular  or,  equivalently,  the  late  of  energy  flux  into  the  crack 
tip  is  nonzero.  flic  criterion  adopted  f  -r  u  tamed  cleav crack  growth  is  that 
the  value  of  this  energy  flux  must  have  a  certain  material  specific  value.  five 

problem,  then,  is  to  determine  the  influence  of  the  intervening  plastic  I  low  in 
screening  the  crack  tip  region  from  the  icm-a/ly  ap;  lied  f  .  i  d  1 1 . .  I  he  analysis  led 
to  conditions  that  must  be  sati.ilTd  foi  .u  lam.  I  teniv  clc.iv.r.c  ciacl.  mouth  in  the 
materials  of  interest.  The  conditions  depend  on  the  I  low  ehaiacterisiics  ol  the 


material,  the  crack  tip  .peed,  and  the  level  of  applied  loadin'/. 


Of  particular 


interest  is  the  result  that,  for  given  material  and  temperature,  there  is  a  minimum 
crack  driving  force  below  which  the  propagation  of  a  sharp  crack  cannot  he 
sustained  under  any  conditions;  and  this  minimum  toughness  is  interpreted  as  the 
crack  arrest  toughness  of  the  material  |MATAGA,  FREUND,  and  HUTCHINSON, 
1987], 

The  crack  propagation  model  has  also  been  examined  from  the  dislocation 
dynamics  point  of  view.  The  interaction  of  a  rapidly  propagating  mode  I  crack  with 
pre-existing  dislocations  is  studied.  First,  the  interaction  of  a  steadily  propagating 
crack  with  a  single  edge  dislocation  near  the  fracture  path  is  studied.  The  stress 
field  of  the  crack  induces  a  shear  stress  on  the  glide  plane  of  the  dislocation  and,  if 
this  stress  is  sufficiently  large,  the  dislocation  will  move  on  its  glide  plane  during 
the  interaction.  On  the  basis  of  a  linear  viscous  drag  model  for  the  motion  of  the 
dislocation,  the  energy  dissipated  in  this  interaction  is  computed.  The  total  energy 
dissipation  per  unit  crack  advance  is  estimated  for  a  distribution  of  pre-existing 
dislocations  by  superposition.  The  end  result  is  an  energy  balance  that  is  similar  in 
form  to  that  obtained  by  the  continuum  viscoplasticity  analysis.  This  new  result, 
although  not  exact,  permits  a  re-interpretation  of  the  continuum  results  in  terms  of 
dislocation  dynamics  concepts. 

The  foregoing  analysis  was  directed  at  understanding  fracture  of  those 
materials  that  may  exhibit  strain  rate  induced  cleavage  but  that  fail  in  a  ductile 
manner  under  slow  loading  conditions.  Of  course,  strain  rate  sensitivity  of  the 
material  response  can  also  influence  the  fracture  response  during  crack  advance 
controlled  by  a  ductile  separation  mechanism.  This  matter  was  addressed  in  a 
related  study  in  which  the  problem  of  steads  growth  of  an  antiplane  shear  crack  in 
a  strain  rate  sensitive  clastic-plastic  material  was  considered  [YANG  and  FREUND, 
1986].  It  was  assumed  that  the  material  is  cl.r.tic-perfcctly  pi. tic  under  slow  loading 


conditions,  but  that  the  inelastic  strain  i.ttc  i,  pi upon iona  1  to  the  difference  between 
the  stress  and  the  y  ickl  tic  tai  cd  to  ..one  power.  In  earlier  w..rk  i.n  this 
problem  reported  b>  "titers  the  po^ ~ i h i 1 1 1 \  of  an  clastic  region  in  stress  'pace  was  not 
considered.  They  concluded  that  the  asymptotic  crack  tip  field  is  completely 

autonomous;  however,  the  local  field  could  not  be  shown  to  reduce  to  a  generally 

accepted  rate  independent  limit  as  the  rate  sensitivity  of  the  material  vanished.  It 

has  been  shown  in  our  work  that  if  the  possibility  of  clastic  unloading  is  admitted 
in  the  formulation  then  the  asymptotic  crack  tip  solution  does  indeed  approach  the 

correct  rate  independent  limit  as  the  rate  sensitivity  of  the  material  vanishes. 

Furthermore,  the  crack  tip  field  loses  the  feature  of  complete  autonomy  under  these 

conditions.  That  is,  the  crack  tip  field  involves  a  free  parameter  that  can  be 

determined  only  from  the  remote  fields. 

2.3  Plasticity  at  High  Strain  Rates 

2.3.1  Pressure-Shear  Plate  Impact 

The  high-strain  rate  pressure-shear  plate  impact  experiment  was  developed  at 
Brown  under  previous  ARO  support.  Under  the  present  grant  this  experiment  has 

been  extended  to  shear  strain  rates,  greater  than  l()',s*1.  which  are  approximately  an 
order  of  magnitude  larger  than  those  obtained  previously.  This  increase  in  strain 
rate  has  been  achieved  by  developing  a  lapping  procedure  that  allows  specimens  to  be 
lapped  flat  and  parallel  to  thicknesses  as  mall  as  35j/m  (KL.OPP.  19,30].  Using  this 
technique  S.  Huang  and  R.J.  Clifton  have  conducted  a  relatively  thorough 
investigation  of  the  strain  rate  sensitivity  of  Oil  1C  copper  up  to  .train  i.ttes  of 
lOV1.  The  shear  stress  increases  with  .train  rate  over  the  entue  tango  ol  train 
rates.  As  one  means  f.-r  under, landing  this  increase,  we  arc  -lev  eloping  the 
capability  to  change  the  strain  rate  during,  the  pressure-shear  experiment.  Currently, 
we  are  trying  to  reduce  the  -.train  i  ate  alter  approximately  O.sjrs  by  using  a  layered 


fiver  plate  consisting  of  a  steel  I  ront  plate  backed  by  a  high-strength  aluminum  plate. 
The  low  impedance  aluminum  plate  causes  reflected  waves  to  partially  unload  the 
specimen  and  thereby  reduce  the  strain  rate  in  the  specimen.  Such  experiments  arc 
expected  to  provide  important  new  insight  into  the  relative  roles  of  the  strain  rate 
history  and  the  current  strain  rate  in  causing  the  strong  increase  in  flow  stress  with 
increasing  strain  rate  at  strain  rates  greater  than  lO's"1. 
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